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. Abstract: Minimal string theory has a number of FZZT brane boundary states; one 

for each Cardy state of the minimal model. It was conjectured by Seiberg and Shih that 
all branes in a minimal string theory could be expressed as a linear combination of the 
brane associated to the identity operator of the minimal model with complex shifts in 
the boundary cosmological constant. Subsequently it was found that this identification of 
FZZT branes does not hold exactly for some cylinder amplitudes but was spoiled by terms 
that are associated with vanishing worldsheet area and are therefore non-universal. In 
this paper we investigate this claim systematically, using both Liouville and matrix model 
methods, beyond the planar limit. We find that the aforementioned identification of FZZT 
branes is spoiled by terms that do not admit an interpretation as non-universal terms. 
Furthermore, the spoiling terms as computed using the matrix model are found to be in 
agreement with those coming from Liouville theory, which also suggests that these terms 
have universal meaning. Finally, we also investigate the identification of FZZT branes by 
replacing the boundary state with a sum of local operators. We find in this case that the 
brane associated with the identity operator appears to be special as it is the only one to 
correctly reproduce the correlation numbers for bulk operators on the torus. 
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1. Introduction 

Minimal string theories are obtained by interpreting Liouville theory coupled to a (p, q) 
minimal model as the worldsheet theory for a string. Advances in the last decade § ||] Q 
have identified the consistent boundary states of the Liouville theory and hence, in the 
string interpretation, allowed the branes to be constructed. The boundary conditions in 
Liouville theory can be classified into two families, FZZT and ZZ. The FZZT boundary 
states |ct)fzzt are labelled by a continuous parameter, a, which is related to the boundary 
cosmological constant. The ZZ boundary conditions, which we will not be considering in 
this paper, are labelled by a pair of integers (m,n). To obtain a brane in minimal string 
theory we must also specify a boundary condition for the matter theory. The consistent 
boundary states \k,l) for the minimal models are in one-to-one correspondence with the 
primary fields and so are labelled by the Kac indices (k,l). To obtain a brane in the 
minimal string we simply take a tensor product of these two states and so, for example, 
the FZZT branes are given by |cr; k, I) = |ct)fzzt ® \k,l). 

Once amplitudes for a disk with FZZT boundary conditions and one vertex operator 
in the bulk had been computed it was noticed in pi that the construction just outlined 
appears to over count the number of FZZT branes. Up to null states it seems that 

W;k,l)= £ £ \ a + i mpl iltl) . (i.i) 

\/PQ 

n=-(H),2m=-(i-l),2 v 
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In effect the disc amplitudes behave as though the information about the matter sector 
boundary condition can be subsumed in complex shifts of the boundary cosmological con- 
stant which describes the Liouville sector boundary condition. It was noted in Q that the 
deviation from the Seiberg-Shih relations for some cylinder amplitudes is non-zero but it 
was argued that the deviation should be discarded as it only appears for certain disk areas 
or boundary lengths. 

In this paper we test the conjecture fll.ip systematically, first using general cylinder 
amplitudes where we show that the prescription described above for discarding deviations 



is inadequate; this case is discussed in greater detail in a recent paper |24j| . We then 
conjecture a new prescription which is correct at cylinder level and test it by considering 
amplitudes for disc with one handle. Disk amplitudes at higher genus have not been 
calculated in Liouville theory so we turn instead to the matrix model formulation. In 
general this approach suffers from the difficulty of constructing the full family of FZZT 
branes which arise from the FZZT boundary state tensored with any matter boundary 



state, although there has been recent progress in this direction [25, pq| . There is however 
a special case where the construction of all FZZT branes is easy; the (3, 4) minimal string 
which corresponds to the critical point of the Ising model coupled to Liouville gravity. 
We use this to show that the conjecture is in fact false. We also show that the failure of 
our conjecture is equivalent to being unable to replace branes other than the fixed spin 
boundary with a sum of local operators. As a by-product of our investigations we show 
there is an easy way to extract correlation numbers from loop amplitudes and verify its 
correctness by reproducing the results of [|ll]] for the (2, 5) minimal string. 

This paper is arranged as follows. In Section || we review the (p, q) minimal string, 



FZZT brane, and the motivation for the identification (1.1) in addition to introducing our 
conjecture. In Section || we examine the failure of the Seiberg-Shih relations for (p, q) 
cylinder amplitudes with arbitrary FZZT brane boundary states and then motivate our 
conjecture. In Section |] we review the calculation of matrix model amplitudes and the 
results we need for cylinder and disc-with-handle amplitudes. We then examine the Seiberg- 
Shih deviation for the disc-with-handle amplitudes and test our conjecture. In Section || 
we investigate the deviations by expressing boundary states as sums of local operators and 
then studying how the Seiberg-Shih relations act on this expansion. Finally, given that our 
results indicate that the different branes are not equivalent, we discuss in the conclusion 
whether we can identify in the matrix model all the dual branes. Various technical results 
are given in the appendices. 

2. (p, q) Minimal Strings and FZZT branes 

The Liouville theory on a manifold with boundary is defined by the Lagrangian, 

Sl = ^~ [ d 2 x^ (g ab d a ^d b 4> + Q<f>R[g] + 4vr M e 2 ^) + f fi B e^dx (2.1) 

47r JM V ' JdM 

where Q = 6 _1 +6, /i is the bulk cosmological constant and \xb is the boundary cosmological 
constant. The primary fields of the theory are V a = e 2a( ^. The worldsheet theory for the 
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minimal string is obtained by coupling Liouville theory to a (p, q) minimal model; this 
requires that b 2 = ~. 

The physical operators of the (p, q) minimal string were first given by [15 as U n = 
O n V a (n) where, 

p ~t~ q — ti 

a(n) = b n 7^ mod p ji^O mod q n£Z (2-2) 

2p 

and O n is an operator containing contributions from the ghost and matter sector in ad- 
dition to derivatives of the Liouville field. An important physical operator is obtained by 
gravitationally dressing the primary fields O rs of the minimal model; such operators are 
referred to as tachyons and take the form, 

T r ,s= ccO r , s V ar , a (2.3) 

p + q-\qr-ps\ 
«r,s = 7T^= 2.4) 



where c and c are ghost fields. The full conformal bootstrap of Liouville theory on the disc 
was performed in [||||], resulting in the identification of a consistent boundary condition 
for the Liouville theory; the FZZT brane. As part of the conformal bootstrap, the bulk 
one-point function on a disc, i.e. the bulk one-point function with FZZT boundary, was 
calculated: 

^(P) = (v P \a) FZZT = (7r^(b 2 )r lP / b T(l + 2ibP)T{\ + 2iP/b f OS{2 l aP) (2.5) 

where a is defined implicitly by, /Ucosh 2 7r6<7 = /x^ 2 sin(7r& 2 ) and \vp) is the primary state 

associated to Vq/2+iP- I n the minimal string the full FZZT brane is \a\ k,l) = \o)fzzt ® 

\k,l), where \k,l) is the (k,l) Cardy boundary state of the minimal model. We will find 

it convenient to follow in rescaling the bulk and boundary cosmological constant so 

that 7rn r y(b 2 ) — > fi and ^B\/^"y(b 2 ) sin(-7r& 2 ) — > fxp- Furthermore, it was noted in § || 

that the results of the conformal bootstrap on the sphere are invariant under the duality 
i 1 

transformation b — > ^, [i — > jl = fib 1 . This self duality is present in the bootstrap on the 
disc if the boundary cosmological constant transforms as, [is A_B; where ^6 = cosh 
The transformed boundary condition is referred to as the dual FZZT brane and provides 
a physically equivalent description of the FZZT brane. 

The amplitude in the minimal string for the insertion of a tachyon in the bulk of a disc 
is obtained using (2.5) M, 



{%,s\cr;k,l) = A 



D 



, irt\ ( Trt 
sin | — sin ( — 
p ) 



T{2ba-b z )T{2b^a-b- 2 -I) x (2.6) 



irt\ ( 7rA _l / -Kta 
Uk-i cos — Ui-i cos — fi 2 p cosh 



p J V q J \Vpq~ 



where t = qr — ps, 2a = Q 1=, is a constant independent of k, r, s and I and 

U n (x) is the nth Chebyshev polynomial of the second kind. This expressions for the disc 
amplitudes motivated Seiberg and Shih to conjecture in M that the construction just 
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outlined overcounts the number of FZZT branes and that, up to BRST null states, the 
following identification holds, 



k;M >= £ £ k + ! !ra U >. (2.7) 

n=-(fc-l),2m=-(Z-l),2 

This relation, which we will refer to as the SS (for Seiberg-Shih) relation, was originally 
obtained by inspection of ( |2.6| ) but was later derived for discs using the ground ring in 
Q. Essentially it states that there is only a single FZZT brane \a; 1, 1) and that all the 
others are related to it by complex shifts in the boundary cosmological constant which has 
somehow absorbed all the information about the matter boundary condition. 

The validity of the SS relation has been checked for the discs and cylinder amplitudes 
|7f||l4j]|p4} and is known to fail in some cases of the latter. For our purposes it is useful to 



define the deviation from the SS relations for a particular amplitude A(k,l;a\X) as 

fe-i l-l 

AA{k,l;a\X)=A(k,l;a\X)- ^ A(l, 1; a + i mp Z^ q \X) (2.8) 

n=-(fc-l),2m=-(Z-l),2 

where the presence of other boundaries and operator insertions is denoted by X. It was 
first argued in [[jj that the deviations noticed there should be discarded as they have the 
following properties: 



Property 1 (PI): The dependence of the deviation on the bulk cosmological constant 
could be cancelled by setting the IR cutoff (that was needed in order to compute the 
amplitude) to the volume of the Liouville direction. 

Property 2 (P2): The inverse Laplace transform, with respect to all boundary cosmo- 
logical constants, of the deviation minus any regularisation dependent parts is zero almost 
everywhere i.e. it is supported only at points. 



Essentially PI says that it can be arranged for the deviation to occur only for surfaces 
of zero area while P2 says that the deviation is non-zero only for particular values of the 
boundary length; and it is argued that in neither case can the deviation have any physical 
meaning. 

In Section ^ we examine the SS relations for cylinder amplitudes with arbitrary FZZT 
brane boundary states. These amplitudes are given in [jl4|] in a rather unwieldy form; we 
give them in Appendix [B] in a more usable form similar to that given recently by [ 24 ] . We 



will show that there exist choices of boundaries for which the deviation possesses neither 
PI nor P2. The SS relations must therefore be modified or discarded. This leads us to 
define the property: 

Property 3 (P3): The only worldsheet geometries that contribute to the deviation are 
those with particular values of boundary length (i.e. zero) and their duals. 
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P3 means that all deviations may be written as the sum of two terms that are dual 
to one another such that the inverse Laplace transform of each term with respect to the 
boundary cosmological constant and dual boundary cosmological constant respectively, 
vanishes almost everywhere. We conjecture that all deviations possessing P3 should be 
discarded (or equivalently that the SS relations for amplitudes hold up to deviations having 
P3). In the rest of this paper we investigate this conjecture. 

3. FZZT-FZZT Cylinder Amplitudes 

In this section we use the cylinder amplitude between arbitrary FZZT branes, which we 
denote Z(k,l;ai\r,s;a 2 ) to show that neither PI nor P2 hold for all cylinder amplitudes. 
First consider a specific example, the case of the (3, 4) minimal string. The general ampli- 
tude in Appendix [B] gives the following results 1 , 

Z(l, 1; ^11, 1; a 2 ) = Z(2, 1; a x \2, 1; a 2 ) = ln(^^) + —L- , (3.1) 

xi - x 2 2V3e 

Z(2,1;<7i|1,1;ct 2 ) = -Info + z 2 ) + ' 



4V3Y 

Z(l,l;oi|l,2;<7 2 ) = Z(2, 1; a x \l, 2; a 2 ) = - ln(-l + 2z\ + 2\p2 Zl z 2 + 2z\) + ' 



2^/3e 



Z(l,2;a 1 |l,2;a 2 )=ln( / ■ Z \,* 2 -) + 



'fo + z 2 )(xi - X 2 ) Ay/^e 

where e is a regulator introduced to cut-off an IR divergence present in the integral and we 
have introduced the following notation which will be useful in the remainder of the paper, 

z = cosh ( ) , x = —^L = cosh ( ^ ) , x = ^= = cosh ( — — ) . (3.2) 

\VmJ \\fpqj \VpqJ 

Clearly we need to remove the dependence on the regulator and in two possibilities 
were considered; either throw away the terms that diverge as e goes to zero (this gives 
agreement with the matrix model), or set it equal to the volume of the Liouville direction. 
If we choose to do the latter then PI should apply, if we choose the former then P2 should 
apply. As an example of a deviation which satisfies PI and P2, consider, 

/3 

AZ(2,l;ai\l,l;a 2 ) = -X_ - l og (xi + x 2 ). (3.3) 

Choosing e equal to the volume of the Liouville direction, \ = \ log where A is a constant 
and recalling that x = (j,b/ we find 

A^(2,l; f 7i|l,l; ( T 2 ) = -- log A - log(/x Bl + // B2 ). (3.4) 



1 These expressions are correct up to an unimportant common normalisation constant and additions of 
numerical constants. 
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Since this is independent of [i it was interpreted as non-universal and therefore set to zero. 
To see that fl3.3Q also satisfies property P2, we take the inverse Laplace transform to get 

_^( Ll )5(L 2 ) + S{In - £ 2 )^-((7 + logii)0(Li)), (3.5) 



AZ(2,1;L 1 |1,1;L 2 ) 



where 7 is the Euler-Mascheroni constant and 9{t) the step function. 

Before moving on to examples for which neither PI are P2 true it is important to 
note that the sum in the Seiberg-Shih relation does not respect the symmetry of the Kac 
table and so for a given amplitude there are two possible Seiberg-Shih relations that might 
possess PI or P2. However, there are cases where both possibilities lead to deviations that 
do not possess either property PI or P2. For an example of a deviation that does not 
possess property PI consider 



A2(l,2;<7i|l,l;<72) 



In 



-! + 2xi 2 + 2\/2a;iX2 + 2x2 2 



2{xi + x 2 ) 



or, since Z(2, 2; ai\l, 1; 02) = Z(l, 2; ai\l, 1; 02), 



AZ(2,2;<7i|l,l;a 2 ) 



In 



-1 + 2xi 2 - 2V2x lX2 + 2x 2 2 ) 



(3.6) 



(3.7) 



where we have suppressed the e dependent term. It is clear that since the argument of 
both expressions is not a homogeneous polynomial in x the dependence on \i cannot be 
factored out and placed in a separate term. Hence the IR cutoff cannot be chosen to cancel 
all dependency on the bulk cosmological constant. For an example of a deviation that does 
not possess property P2 consider, 



AZ(l,2;<7i|l,2;<7 2 ) = hi 



X\ - X2 



2(xi - x 2 ) 



or, since Z(2, 2; ai\2, 2; cr 2 ) = Z(l, 2; o\ |1, 2; cr 2 ), 

(xi - x 2 ){x\ + x\ - 1) 



A£(2,2;ai|2,2;a 2 ) = In 

;e Laplace ti 
expression for the amplitudes, 



(3.8) 



(3.9) 



2(xi - x 2 ) 

To take the inverse Laplace transform of these deviations it is easiest to use the integral 



Z(k,l;ai\r, s;a 2 ) 



2tt 2 
V2pq 



dP cos(2vrc7iP) cos(27rcr 2 P) sinh(^) _ 2 ~/ P 



P 



sinh(2vr6P) sinh(^) 



pq 



where, 

A p (fe,r) A 9 (i,s) p— 1 (j— 1 

j ? * A r,.w= E £ E E 

jj=|fe-r|+l,2 p=|i-s|+l,2 a=l b=-(g-l) 



. 7Tt . 7Tt 

sin — J sm( — ) x 
P Q 



(3.10) 



(3.11) 



irt s 



C/ ?7 _i(cos — )U -i(cos — 
' p F q ■ cos 



7rt 
P9 



COS z 
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and \ p (k, r) = min(A; + r — 1, 2p — 1 — k — r). This integral is a generalisation of the integral 
appearing in |23| and [[7j; its derivation may be found in Appendix || [24]. Following [23] 
the inverse Laplace transform may be computed by noting that 



7r6cos(27rPcr) 



/ ^e- MUosh ^ bt7) K2 1 p{Ml) (3.12) 
Jo I b 



2Psinh(^) 

and then substituting this relation into ( |3.10|) . The result is 

/•oo sinh(^£) 2-niP 

C^W^l'Mr^-M)*- / dPtpp(Lx)il)p(L 2 ) , ^ F ( ) (3 . 13) 

Jo smh(27rP6) ^/pq 



where ipp(L) = \J P sinh(2-7rP '/b)K2iP/b(L). By noting that the integrand is symmetric, the 
integral may be computed by closing the integral along the entire real axis by a semicircular 
contour. Since the Seiberg-Shih deviation is merely a linear combination of amplitudes its 
double inverse Laplace transform with respect to each boundary cosmological constant also 
exists and is given by 

roo s inh(^£) 2m P 

£- 1 [AZ](A:,;;L 1 |r, S ;L 2 )« / dP^ P {L{)^p{L 2 ) . ^ AF M , r , 8 (-— ) (3.14) 

J smh(27rPb) ^/pq 

where 

AF k; i ir , s (z) = F k j^ s (z) - U k -i(cos qz)Ui-i(cos pz)Fi^ s (z) (3.15) 

and we have assumed that the Seiberg-Shih transformation has been applied to the (k, I) 
boundary. An important property of Ai^ j r s (^) is that it is an entire function and therefore 
the only contributions to the integral comes from the poles at P = ^ where n£Z. It is 
now easy to show PI and P2 are not satisfied for (|3.8| ) and ( |3.9| ); in both cases there exist 
poles when P = |r and therefore when the integral contour is closed around them it will 
result in a sum of terms containing Bessel functions of the form Kin (L) which have global 

3 

support in L. 

We now want to motivate the idea that P3 is a reasonable extension of the criteria for 
throwing out terms that spoil the Seiberg-Shih relations (and indeed is the most generous 
we can think of). Although for the moment we use the cylinder amplitudes as a motivation 
for introducing P3 we will later analyse the case of the disc-with-handle amplitude more 
carefully. 



Consider the integral representation of the amplitude (|3.10[) ; it may be computed by 
extending the region of integration to the entire real line and then splitting the integrand 
up into terms for which the contour may be closed in either the upper or lower half plane. 



Assuming that o\ > then upon substituting in ( |B.15| ) and (B.16) this results in the 
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expression, 

Z(k,l;a 1 \r, s;a 2 ) oc F k i r JQ)(^- - 71-01)— != + (3.16) 

2e jpq 

EApq . to to to , to 7rn0- 2 , 
(7fc_i(cos — j(7 r _i(cos — jtV;_i(cos — jt7 s _i(cos — )e v?? cosh h 
n p p q q Upq 

71 = 1 ' v 

nj^Omodp 
n^Omodq 

El, . _ ^ .to, sin(7rn/g) 
n q sm{irpn/q) 

n^O mod ij 

E -(-iriW-)^^e^ cosh(TO, 2 /6) + 

n^O mod p 

E^(-ir^F to (TO)e-^cosh^. 
^— f ring 

The above expression may be understood as representing a cylinder amplitude as a sum 
over disc amplitudes with various local or boundary operators inserted. If we compute 
the deviation using this expression the first sum cancels while the others remain with the 
replacement Fk t i <r ,s ~^ A.Fk t r s . The sum in which each term is proportional to cosh (-7T 67102) 
arises from the poles in (3.1C) due to the factor of sinh(^p) in the denominator; it is 
precisely these terms which vanish 2 when we take the double inverse Laplace transform 
with respect to the boundary cosmological constant. These terms have the interpretation 
of some form of descendant dual boundary length operator inserted on the boundary of a 
disc. On the other hand, the sum in which each term is proportional to cosh(TO<r 2 /6) arises 
from the poles in ( |3.10| ) due to the factor of sinh(27r6P) in the denominator; it is these terms 
which, under a double inverse Laplace transform with respect to the boundary cosmological 
constant, become the terms proportional to Bessel functions of the form K±n (£2). Similarly, 
if we were instead to take the double inverse Laplace transform with respect to the dual 
boundary cosmological constant, the terms proportional to cosh(7rncr2/6) would vanish 
and the terms proportional to cosh^nfro^) would become terms proportional to Kin (L 2 ). 
We therefore see that the deviation computed from ( 3.16| ) can be understood as a sum of 
terms which vanish under a double inverse Laplace transform with respect to the boundary 
cosmological constant or its dual. 



4. Testing the Seiberg-Shih relations using a Matrix Model 

In the last section we argued that the deviations for all cylinder amplitudes have property 
P3. Obviously it would be interesting to know whether this is true for deviations of any 
amplitude. Unfortunately amplitudes more complicated than the disc and cylinder have 
never been computed using the Liouville approach. An alternative is to use the matrix 
model formulation of minimal string theory in which the computation of amplitudes with 
arbitrary numbers of boundaries and handles is straight-forward. 
2 Actually they have point-like support if we were being more careful. 
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The partition function of a hermitian iV-matrix model can be interpreted as generating 
graphs with vertices carrying iV labels with weighting determined by the action of the 
model. By tuning the parameters that appear in the action we may reach a second order 
phase transition which will be described by a conformal field theory. By studying the 
critical exponents, the CFT associated with the phase transition can be identified with 
Liouville gravity coupled to a minimal model. The disadvantage of using the matrix model 
is that it gives the amplitude in a fully integrated and summed form and so it is hard to 
understand the structure of the amplitude in terms of continuum concepts such as states 
circulating in a loop. Furthermore, generally it is obscure how the graph labels map to the 
conformal field theory degrees of freedom and even more importantly, the matrix model 
appears to not contain all the boundary states present in the minimal string. Indeed, this 
was used as evidence in Q in support of the Seiberg-Shih relation. However, there are 
special cases where the relation is manifest; the (3, 4) minimal string admits a formulation 
as a matrix model in which it is easy to construct all boundary states. This is the model 
we will study in this section. In recent work of [p^, 26] a number of other boundary states 



of the (p, q) minimal string have been constructed in the matrix model formulation and it 
would be interesting to extend our results to those cases. 

The (3, 4) minimal model coupled to Liouville gravity may be thought of as the CFT 
that describes the critical point of an Ising model on a random lattice. It is easy to 
formulate a 2-matrix model in which the weighting of the graphs allows the two labels to 
be interpreted as the spin degrees of freedom in the Ising model with a coupling dependent 
on the coupling in the matrix model action. Such a matrix model is 

Z = J dM x dM 2 exp (-j Tr (-cM^ + ^(Mf + M|) + i(Af? + M 2 3 )^ ^ (4.1) 

The Feynman rules for this theory generate graphs of coordination number three with each 
vertex generated by either ^Mf or |-Mf . Such a graph may be interpreted as a discrete 
surface composed of triangles with each vertex at the centre of each triangle and with label 
spin + or — depending on if the vertex is generated by Mi or Mi. The coupling between 
neighbouring spins is controlled by the parameter c whereas g controls the cost associated 
with adding more vertices to the graph. It is g and c we use to tune the matrix model to 
its critical point which for our matrix model we do by setting, 

9 = ffc(l - a 2 w), c = c c , (4.2) 

where r\ is some constant which can be determined by comparing with the Liouville theory 
and the critical point is achieved by letting o — > 3 . We will often refer to this limit at the 
scaling or continuum limit. 

The conformally invariant boundary conditions of minimal CFTs are in one-to-one 
correspondence with the primary fields of the theory. For a CFT that describes the critical 
point of some discrete model the boundary conditions can sometimes be understood as 
universality classes of boundary conditions in the discrete model. As we approach the 



3 The constants g c and c c have the well known values g c = 10c;? and c c = -L(2\/7 — 1). 
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critical point of a discrete model all boundary conditions in a given class will flow to the 
same boundary condition in the CFT. In particular for the (3, 4) minimal string its three 
conformal boundary conditions are the continuum limit of the discrete configurations of 
boundary spins consisting of all spin +, all spin — , or free spins. It is these boundary 
conditions we want to implement in the matrix model. 

The usual way of inserting a boundary in the matrix model is to compute the resolvent, 

WMl (x) = /itt-J-N S 1 V (4.3) 



N x-Mi N ^ x r - 

The graphs that contribute to each quantity (TrM n ) correspond to triangulated surfaces 
with a boundary of length n composed of vertices generated by M\. The function Wm x is 
a generating function for such quantities. We will adopt the following notation for more 
general quantities, 

W Ql Fl Q2 F 2 .., Q3 Gl Q4 G 2 ..,...(xi, x 2 , y 1: y 2 , ...;...) = 
1 / Qi(M 1 ,M 2 ) Q 2 (M l ,M 2 ) \ 

N Axi -F 1 (M u M 2 )x 2 -F 2 (M 1 ,M 2 ) "') "' 1 ' ' 

1 T (, Q 3 (M 1 ,M 2 ) Qa{M u M 2 ) 
- lr ' ' 



N y yi -G l (M 1 ,M 2 )y 2 -G 2 (M l ,M 2 ) 



By tuning the . . . as we take the matrix model to its critical point we can extract 

continuum quantities corresponding to amplitudes with macroscopic boundaries. By way 



of example, for Wmx(x) we use ( [4.2[) and set, 

x = x c (l - a d Kfi B ) (4.5) 

where d is chosen to produce a non-trivial limit and k is chosen to agree with Liouville 
theory. As we let a — > 0, Wm% will have an expansion in powers of a 4 , 

^(1,(1 - a d ¥B )) = ^a d 'f,( Wl/ i) + ^%> B ,ri I b.t . (4.6) 

i 

In this expression Wi(^B, i s analytic in both arguments whereas Wm 1 is defined as being 
the first term non-analytic in /i# and /x. It is Wm\ that corresponds to the continuum 
quantity and so in this case is the partition function of the continuum theory defined on 
a surface with a finite sized boundary with the boundary condition determined by which 
universality class the discrete quantities (TrM") belong to. This is slightly complicated by 
the fact that the resolvent actually corresponds to a boundary with a marked point and 
therefore must be integrated with respect to the boundary cosmological constant before 
comparing to Liouville theory. This motivates introducing the integrated quantity, 



4 When taking the scaling limit in the remainder of the paper we will choose 77 = 1 and k = yj g c /c c . It 
should be kept in mind though that by changing r\ and k we can renormalise /is and /t. 
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where we have allowed for more general amplitudes which have more than one boundary. 



The generalisation to amplitudes of the form (4.4) is obvious. 



For the matrix model (4T) it is easy to construct resolvents which generate boundary 
conditions that flow to the three different boundary conditions in the CFT: 



W Ml (x) 
W M2 (x) 

W Ml +M 2 (x) 



\ Tr - 1 



N x- Mi 

1 m 1 
:Tr- 



N x - M 2 



(4.8) 



1 



:Tr 



1 



N x - (Mi + M 2 ) 



The first two resolvents Wjv;, generate surfaces on which only Mj vertices appear on the 
boundary. Since the Mj vertices map directly to the spin degrees of freedom these flow to 
the fixed spin boundary conditions in the continuum limit. The final resolvent generates 
boundaries in which both types of vertices appear with equal weighting. This will therefore 
flow to the free spin boundary condition. We can parameterise these resolvents in the 
following way, 



X 



X 



a 



Mi 



[l + a) 



Mo 



(4.9) 



This generates graphs with weighting corresponding to the Ising model on a random lattice 
but with a boundary magnetic field, to which the parameter a is related. Clearly, we can 
choose a = ±1 to reproduce the resolvents Wmi- By taking a to infinity we may also 
obtain the resolvent Wmi+M^- One method to compute the resolvent Wx is to make the 



following change of variable in (4.1), 



S = Mi + M 2 + 1 + c, 



X 



(1 



a) 



-Mi 



Mo 



The partition function in these new variables takes the form, 



where, 



Z = / [dSdX] exp 



1 



N 
9 



Tr (X 2 S + aXS 2 + V(S)) 



(4.10) 



(4.11) 



V(S) = -(l + 3a 2 )S 3 



C ^ 2 1 



-S 2 + -(3c 
2 4 V 



l)(l + c)5. (4.12) 
Generic values of a (i.e. the Ising model on a random lattice in the presence of a boundary 



magnetic field) were investigated using matrix model techniques in [18], in which the disc 
amplitudes were calculated. The method employed in [|ll| to solve the matrix model is 
combinatorial and so is not easily generalised to compute more complicated amplitudes 
and so we will not employ it here. However, it is clear that if we are only interested in 
amplitudes in which all the boundary conditions are of the free spin form we may set a = 



in ( 4.11 ) to obtain the 0(1) model whose solution is well known. Similarly if we are only 
interested in amplitudes in which all the boundary conditions are of the fixed spin form 
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we may simply use ([O]). This is in fact the case for most of the results we require for 
the remainder of this section and so they can be obtained form the literature 0][17]. For 



more complicated amplitudes involving both fixed and free spin boundary conditions it is 
necessary to solve the matrix model ( [4.1 1| ) ; this may be done using loop equations and we 
give details of this in Appendix |Cp. Although we may look up many of the results in the 
following section in the literature, we will, in order to discuss all the matrix model results 
in a common language, phrase the discussion as if we had computed the various resolvents 
from ( 4.11| ) using the loop equations. 



To take the scaling limit of Wx at a = —1 we may use the results of || or the results 
in Appendix [C]. The critical value of x is x = — c and the solution to the loop equation 
gives, in the scaling limit, 




For we may use [[jj]] or solve ( |C.16 ). The critical value of x is at x = and computing 

the scaling limit of we get, 




Wm 1+ m^b^) = --^t- \[^ + \ ^-M + Fi-V^-H • (4-14) 



Here we run into one final technicality. For a resolvent that produces graphs with some 
configuration of spins on the boundary Wf(Mi,m 2 ) > there exist an entire family of resolvents 
which produce the same configuration of boundary spins but with different overall weighting 

W ±F(M U M 2 )(. X ) = XW F(M U M 2 )(^), (4-15) 

where A is a constant. Changing A has the affect of renormalising //g, however we want the 
same normalisation for all continuum quantities and so this fixes A to be v2 for Wm 1 +m 2 - 
Therefore using ( [4.14 ) and (4.13), we have that 



UJ (M 1 +M 2 )/V2^B,fJ') = V2u Ml (fJ, B ,}J,). (4.16) 

This is precisely the relation between these two disc amplitudes found in the continuum 
calculation. 

One may wonder what the scaling limit of Ifj' ' is for generic values of a; given that 
there are only a finite number of Cardy states the expectation is that there should be no 
new non-trivial scaling limits. For generic a the critical values of x occur at -^r and -%r- 
Taking the scaling limit, for both critical values of x, gives 

W-X/a^B,^) = {flB, AO (4.17) 

which confirms our expectation. We have therefore demonstrated that we can build all the 
boundary states of the (3, 4) minimal string. This will be confirmed further by reproducing 



We have, as a cross-check, computed all the amplitudes in this paper using the loop equation method 
detailed in the appendix. 
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the cylinder amplitudes shortly. The nature of the scaling limit when a = will be 
addressed in the Conclusion (^). 

The one case in which we are interested which cannot be found in the literature is 
Ws-x- This may be obtained by the method of loop equations and we give the details in 
Appendix [C]. The computation results in 

us;x(HBi,liB2,n) = -r^\og{-l + 2zl + 2V2 Zl z 2 + 2zl) (4.18) 

for the scaling form of Ws-x, which is in exact agreement with the Liouville result (up to 
renormalisations of \xb and \x). Of course the other cylinder amplitudes may computed in 
this matrix model or using B |l7| and we find the scaling forms to be, 



1 Z\ — Z2 

UX;X{VBl,HB2,M = Y0^| x 1 -X 2 ' 

WX ; y(MB1,MB 2 ,M) = ~J^2 lUZl + Z2 ' ( 4 ' 2 °) 



ws-,s(j*Bi, Mb 2; M) = TFT2 ln 7 1 Z \r ~ v ( 421 ) 

IOC 2 [Zl+ Z2){Xi - x 2 ) 

again in agreement with the Liouville calculations. The fact that these amplitudes agree 
with the Liouville calculations including terms that were classified in Q as non-universal 
suggests they really should have a physical meaning (and shouldn't be thrown away). This 
issue was raised in |]lT[ . 

In order to test if P3 is a property of all deviations we need to examine more com- 
plicated amplitudes. The 1/N correction to the disc amplitude is the amplitude for a 
disc-with-handle. These computations have already been done in [^] and ||, the result 
being, 



(D, m, 1 z(7 - 24z 2 + 48z 4 ) _i3 
u4/>B^) = u^{hb,h) = ~ 2 i 5 | 36c ( _ 1 + 4z 2)3 » 6 ' ( L22i 



(l) v V2 z(5 - 24z 2 + 48z 4 ) _js ft N 

V 1+ M 2 )/^^^) = - 2 | 5 | 36Cc (-1 + 4^)3 6 • (423) 

If we now compute the deviations we find, 

Au;«(mb,m) = 0, (4.24) 

(i) 2^4 z(-3 + 41z 2 - 208z 4 + 416z 6 - 512z 8 + 256z 10 ) _7 

Aw (m 1+ m 2 )/^2^ b '^ - (-1 + 2z) 3 (l + 2z) 3 (l - 16z 2 + 16z 4 ) 3 M 6 ' 

We now want to consider if Aw^ ,-(//#, u) possesses P3. First note that if we 

write 

Aj m + M 2)/ V2^^ = A &)+B(x) (4.25) 
then in order for A and B to have an interpretation as the contribution of certain geometries 
to the amplitude then both A and B must be positive for physical values of their arguments, 
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in particular as x, x — > oo 6 . Furthermore, since 



-11 
2— 1 



A -(Mi + m 2) /VI^^) = - iff* ~ -V4 ( 4 - 26 ) 

this means that A(x) = ax~ 7 ^ 3 + o(x _7//3 ) or B(x) = bx~ 7 ^ + o(x~ 7 ^) otherwise for 
sufficiently large x either A or B would be negative. We therefore conclude both the inverse 
Laplace transforms of A and B with respect to either x or x exist. For Au; , , , , to 

possess P3 we require that .C^ 1 ^] and 7 have point-like support. However it is a 

well known theorem that a function with only point-like support may be expressed as a 
sum of derivatives of ^-functions, 

oo oo 

/(*) = EE a ^ 6(n) (* - l ^m) (4.27) 

771=1 71=1 

and such a function will have a Laplace transform of 

oo oo 

c [/](*) = EE a «.™ sVS '" ,ra ( 4 - 28 ) 

711=1 71=1 

Given that we require that £~ X L4] and /^[.B] are point-like supported this implies A(x) 



and B{x) have an asymptotic behaviour as x — > oo inconsistent with ( 4.24 ) and hence the 
deviation can not possess P3. 

5. The Seiberg-Shih Relation and Local States 

Suppose that the boundaries considered in the previous sections could be expressed as a 
sum over local operators in the theory, such as, 

oo 

(fc,Z;<ri| ~ y~]f n (k,l,(Ti)(On\ , (5.1) 

71=1 

where O n represents some generic local operator and ~ means equivalence up to terms that 
diverge as o\ — > oo. If this were true then the effects of the Seiberg-Shih transformations 
on a boundary would reduce to studying their effect on the above series. The expansion 
of a boundary in terms of local operators was investigated in [ |10| | in which it was used to 
compute the one and two point correlation numbers of tachyon operators on the sphere for 
the (2, 2p + 1) minimal strings. Later it was shown |l2j that such a technique could also 
compute the three point function correctly for the (p,p + 1) models. 



One technical point raised in [10] and [12] was that the boundary state could only 



be expressed as a sum over local operators if we included operators not in the BRST 



This may seem to contradict the expressions given earlier for cylinder amplitudes e.g (3.8) however in 
this case some of the contribution from the singular geometries has cancelled between A and B; if they are 
reintroduced then the resulting A' and B' will be positive. 

7 We denote in the subscript the variable the inverse Laplace transform is taken with respect to. 
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cohomology of [15]. That the operators in the cohomology of [15], which we refer to as 
LZ-operators, are not all the physical operators of the minimal string is easily verified 
by noting that the cohomology does not contain the boundary length operator § dM dxe^ . 
Such an operator can be inserted on the disc by differentiating with respect to the boundary 
cosmological constant, giving the result 

j> dxe b ^ = A(p,q)U k ^ (cos ^ (cos ^ cosh » ( 5 - 2 ) 

where A(p, q) is a constant of proportionality that only depends on p and q. Indeed there 
exists a whole family of such operators with Liouville charge b p+c ^ p nq where n £ Z + which 
we will denote by U nq . 

In order to compute the coefficients f n (k, I, o~\) appearing in ( |5.1[ ) we need to know the 
amplitude for any local operator inserted in a disc with a (k, I) FZZT boundary condition. 
For the tachyon operators the amplitude is fl2.6|) . For the non-tachyon operators we have 
not computed the exact amplitude using Liouville theory, however it is easy to find the 
contribution to the amplitude from the Liouville sector in the semi-classical approximation 
by solving the Wheeler-DeWitt equation, combined with the matter sector we get 8 



(U n „\ k, I; a) oc cosh ( — Jzr ) • (5-3) 

V VM ) 



Finally, if we apply the Liouville duality transformation to the above expression we get 

(7TTipO~ \ 
VP1 J 



where we have introduced the operators U np which have Liouville charge 62±2_™£, The 
expansion (|5.1| ) can now be given a more concrete form, 

oo 

(M;<ti| ~ y"]/n(M)Pi)(Z4| • (5-5) 

n=l 

Having introduced the necessary technical results we now want to prove the following 
claim, 

If all the FZZT branes in the (p, q) minimal string can be replaced by local operators then 
the deviation from the Seiberg-Shih relations is caused only by the non-LZ operators. 

Consider the cylinder amplitude Z(k, I; o~i\r, s; 02). If we replace one of the boundaries 
by a sum of local operators using (5.5) then we get, 



Z(k,l;a 1 \r, s;a 2 ) = ^ f n (k, I, ai)(U n \ r,s;a 2 ). (5.6) 



n=l 



8 In fact the semi-classical approximation has been shown to give the exact answer in the case of LZ- 
operators inserted in the disc; we are going to assume this is also the case for non-LZ operators. 
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By comparison of fl3.16|) with ( |5.6[) , we find for the coefficient of the LZ operators, 
4pqAc 1 vrn 7rn N 



f n (k,l;ai) = -^-Vr- £4-1 (cos— )C/ z _i(cos— )/x *e = f n (k,l)e ^, (5.7) 
A£>(n) n p g 

where Ac is just a numerical constant and n is not a multiple of p or g. Furthermore, note 
that the LZ operator coefficients satisfy 

/„(M;(7i)= E E /n(M)e vmK ! v*> (5.8) 

o=-(*-l),2 6=-(i-l),2 

This shows that if all boundary states admit an expansion in terms of local operators 
then the Seiberg-Shih relations transform the coefficients of the LZ-operators correctly. 
Hence any deviation from the Seiberg-Shih relations must come from the non-LZ operators. 
However, it is the non-LZ operators one would expect to give deviations compatible with P3 
as they correspond to boundary operators and their duals. One should then be suspicious 
that in fact the boundary states may not be expressed as a sum over local operators. We 
will now see this suspicion is borne out by examination of the cylinder amplitudes. 

If we note that for cylinder amplitudes, Z(k, I; cri\r, s; a 2 ) = Z(k, I; cr 2 \r, s; ai) then 
together with (|5.5| ), this implies, 



00 00 

E 

n=l n=l 



f n (k,l,ai){U n \r,s;a 2 ) = y~]fn(r,s,ai)(U„,\ k,l;a 2 ). (5.9) 



Since {U n \ r, s; a) = An' s ^ g n {&) where g n form a set of linearly independent functions then 
this implies 

f n (k,l,a) = A^h n (a), (5.10) 

where h n is some function. We conclude that if all states admit an expansion of the form 
Q5.5| ) then cylinder amplitudes can be expressed as, 

00 

Z{k,l;ai\r,s-a 2 ) = ^A^At ) K(a 1 )g n (a 2 ). (5.11) 

n=l 

We now want to see if the cr-independent coefficients of each term in ( |3,16| ) has the above 
property. Consider the coefficients appearing in the second sum of ( |3.16| ) , 

vrn sin(7rn/q) 

Fk,l,rA — — ( T\ = ( 5 - 12 

q sm(irpn/q) 



\ p (k,r) \(l,s) 

E 

j/=|fc-r|+l,2 (t>=\l-s\+l,2 



t (-ir^ifooB^Ve^-^)-^, 

IM + i,2 \ 9 J\ PJ 

where X p (k, r) = min(/c + r — 1, 2p—l — k — r) and n/0 mod q. We want to know whether 
this can be written in the form An ' A n ' . That this is not true in general is shown by 
considering the case (p, q) = (4,5), for which we have, 

ttti smjirn/q) 

F 2 ,2,2,l( — )—, T\ = °> v 5 - 13 ) 

q sm{irpn/q) 

im sin(7m/g) 1 

*l,l,2,l( 1-T = 7T> ( 5 - 14 ) 

q sm{irpn/q) 2 
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which is clearly incompatible with the factorisation required ( 5.11 ), This lack of factorisa- 
tion appears to be generic; for instance if we consider the (3, 4) model we find that there is 
not even a subspace of states which admit an expansion in terms of local operators. This 
can be seen by considering the matrix, 



vrn sm(irn/q) 

£k,l,r,s{ )- 



q sm^pn/q) 



( %=fl I fcos^ \ 

x 2(-ir | ( _ ir+ i cos 3^ 

\§cos^p §(-l) n+1 cos^ ±(-l) n (3-4cos 2 ^f ) J 



(5.15) 

where the columns and rows correspond to (k, I) = (1, 1), (2, 1), (1, 2). If there exist A^' 1 ^ 
and An' s ^ such that i 7 fc,i,r,s(^ ; ) = An An' s \ then the determinant of the matrix in ( 5.15| ) 
must be zero which is not the case. 

The same conclusions can be reached by studying the other coefficients associated to 
the non-LZ operators in ( |3.16| ) . This leaves us in an awkward position; there is no subspace 
of states that can be expanded in terms of local operators unless we allow the subspace to 
be one-dimensional (in which case why choose a particular boundary as being the one that 



can be expanded in terms of local operators), yet from the results of JO]] and [12] it is clear 
information about insertions of operators can be extracted from loops. It is worth noting 
that from the examination of cylinder amplitudes in which one of the boundaries is a ZZ 



brane it appears that ZZ branes do not couple to LZ states [24, 27, 

An interesting observation is that all correlation numbers were extracted from loops 
with fixed spin boundary conditions. We have found further evidence that correlation 
numbers can be extracted from the fixed spin boundary by computing the 1-point correla- 
tion numbers on the torus starting form the disc-with-handle amplitude in the case of the 
(2,5) model. We have found exact agreement with the computation done in [O]. A trick 
that makes this computation easy is that all amplitudes for the minimal string are easily 
expressible in terms of z, which is the uniformising parameter of the auxiliary Reimann 
surface (IJ or equivalently spectral curve. If we write these amplitudes in terms of a new 
variable, w, defined by z = \{w + w -1 ) we may easily compute the large w expansion of 
the amplitudes. Since z = cosh -^=, the large w expansion will be an expansion in terms of 



the functions which is exactly what is required to compute the correlation numbers. 
Explicitly, if we have an amplitude (1, 1; a\\ X) then 



oo 

(l,l;ai|X) =^,4 n u;- n = X^ n e"W, (5-16) 



n=l n=l 

but we also have 



(l,l;a 1 |X)=^/ n (l,l)e"W(W n |X) (5.17) 



n=l 

giving, 



(U n \ X) = s An . (5.18) 
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The reproduction of the torus correlation numbers in Appendix [A] is evidence that this 
procedure is correct and that the fixed spin amplitude can be used to compute correlation 
numbers of local operators. Furthermore, we can obtain evidence that the non-fixed spin 
boundary states do not admit an expansion in terms of local operators by applying the 
above procedure to the disc-with-handle amplitudes we computed for the (3, 4) model. We 
have the expansions, 



(1) 
Mi 



W (M 1+ M 2 )A/2 (/iB ' M) 



12 7 
2 3 5 3 cA 6 



12 7 
2 3 5 3 cA 6 



1 



1 



24w 
1 



72w 5 
1 



72w 7 
1 



24w 72w 5 



+ 



+ 



17 



72w 
7 



72w 7 72w 



(5.19) 
(5.20) 



If the expansion in terms of local operators is valid for both these boundary states then 
we would expect the coefficients of each term in the large w expansion to be related by a 
factor of \[2 as can been from ( |5.7[ ); this is clearly not the case. The implication of these 
results is that the fixed spin boundary condition is special and is the only one that admits 
an expansion in terms of local operators. 



6. Conclusions 

In this paper we gave examined the Seiberg-Shih relations ( |1.1| ) for cylinders and discs- 
with-handles. For cylinders we reproduced the result of [0]|j2J] that the identification of 
FZZT branes is naively spoiled by terms arising from the gravitational sector of the theory 
and investigated the claim in Q that all such terms are non-universal. We found that the 
terms spoiling the identification cannot always be associated with degenerate geometries of 
the worldsheet as other geometries which are degenerate in a dual sense are also present. 
This lead us to conjecture that the terms spoiling the FZZT brane identification could 
always be interpreted as arising from degenerate and dual degenerate geometries. We 
checked this conjecture by computing the disc-with-handle amplitudes for free and fixed 
spin boundary states in the (3, 4) minimal string using matrix model techniques. We found 
that the deviation from the Seiberg-Shih relations in this case could not be written in 
the way we conjectured. Given that our conjecture was a very generous interpretation of 
what terms might be unphysical we conclude this is strong evidence that the FZZT brane 
identification conjectured in Q does not hold at all levels of perturbation theory. 

We also considered an alternative approach to testing the Seiberg-Shih relations by 
expanding boundaries in terms of local operators. If such an expansion were possible 
we showed that it would lead to deviations from the Seiberg Shih relations consistent 
with our conjecture. We then gave explicit examples for cylinder amplitudes where a 
local operator expansion of boundary states fails. This lead us to a paradox; how could 



correlation numbers of local operators be extracted from boundary states in [10] ]T2[ if such 
boundary states cannot be expanded as local operators? To investigate this we computed 
the disc-with-handle amplitude in the (2,5) model using matrix model methods and then 
showed that the fixed spin boundary condition yielded an expansion from which correlation 



numbers in agreement with the results of [11] could be extracted. This lead us to conjecture 



- 18 - 



that the fixed spin boundaries are special in that they allow an expansion in terms of local 
operators. 

There are many obvious generalisations of the work done here, the most obvious being 
that it would be interesting to extend the above results to the (p, q) case. Since the Liouville 
techniques are not suited to higher-genus computations such a project would have to be 
tackled using matrix model methods. The main problem with this approach is that it is 
difficult to represent all the states for a given minimal string in a matrix model; in the 
above we chose the Ising model precisely because there was an obvious mapping between 
the spin degrees of freedom and the matrix fields. However this problem may now be 
less serious given recent work [25, 26] in which many boundary states were constructed 
in the matrix model formulation. Another strategy for attacking the same problem may 



be the geometric recursion techniques developed by Eynard et al. in [19] as this provides 
an efficient way of computing many amplitudes for a given matrix model. Such methods 
would also allow a relatively easy computation of the cylinder-with-handle amplitudes; 
these would be interesting to study as one could then check if quantum corrections destroy 



the fact, noted in [24, 27, 28], that FZZT-ZZ cylinder amplitudes are consistent with the 
Seiberg- Shih relations. 

Another generalisation would be to perform a similar investigation in other non-critical 
string models. An obvious choice would be non-critical superstrings as they perhaps have 
more physical relevance in addition to being much better behaved. A second and perhaps 



less obvious choice of string model would be the causal string theories developed in j2fj, 21 



22]. Currently these models have only been solved for a target space of zero dimensions 
and so such a project would require a generalisation of the models to include worldsheet 
matter. However, even in the zero dimensional case such causal string models are better 
behaved and many of the odd features of the usual non-critical string are absent due to the 
lack of baby-universe production; in particular the gravity sector in these theories seems 
to be much weaker. Since the terms spoiling the FZZT brane identification are due to 
the gravity sector (at least for the cylinder amplitudes) might they have an interpretation 
in terms of baby-universe over production and if so might they be absent in these causal 
string models? 

Finally, perhaps the most obvious omission in this paper is a discussion of the nature 
of the dual branes. In ||] the two matrix model description of the (p, q) minimal string 
was considered and it was shown that the continuum limit of the resolvent for one of the 
matrices was the disc amplitude with a fixed spin boundary. The continuum limit of the 
resolvent for the other matrix gave the dual brane amplitude. In our matrix model ( 4.1 1) ) 



it is not obvious if we can construct a correlation function which gives the disc with dual 
brane boundary conditions. In fact it is present in the form, 

Wa(x) = I — Tr — t ) (6.1) 

AK ' \N x- (Mi -Mi) I ' 

which may be obtained from the general X resolvent ([O]) by setting a = 0. The resulting 
loop equation is easily solved as it is quadratic in x 2 , however the scaling limit is obtained 
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using a different scaling of x, 

- i 

x=x c + a3fl B y'g c /c, (6.2) 

which results in 




Wife(Mi-M2)(£B,/i) = j-j- /iB + V^i-^ + Mb-V^I - / 1 > ( 6 - 3 ) 



where A; is a numerical constant. If we identify /xb in this expression with the dual boundary 
cosmological constant then this is precisely the amplitude expected for the dual brane 
obtained from the spin + or spin — boundary state. This explains the a = possibility we 
left unexplained at the end of Section |j. We therefore see that the matrix model is able to 
reproduce all the boundary states found in the (3, 4) minimal string apart from the dual of 
the free spin boundary condition. In Q this was enough as they claimed that there is only 
one brane and hence only one dual brane in the theory. However, in light of our results 
there should exist other dual branes in the theory. Are such states present in the matrix 
model and if so can they be represented as some form of resolvent? This question and the 
others outlined we leave to future work. 
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A. Computation of torus 1-point correlation numbers 

The disc-with-handle amplitude for the (2, 2p + 1) minimal string has been computed 
numerous times using matrix model techniques. For the case of the (2, 5) model it takes 
the form, 

WAf(/iB,/i r- — 3 f 1 4 (A- 1 ) 

2234576 z 

The boundary condition on the disc is the equivalent of our spin up boundary as it is com- 
puted using the resolvent of M where M is the only matrix appearing in the matrix model. 
We will now use this to compute the torus one point correlation numbers, for tachyon 
operators, by replacing the boundary with a sum of local operators as in ( |5.5| ). That this 
computation produces results that match exactly the results of [|ll]] is evidence that we may 
replace this boundary by a sum of local operators. Obviously, in order to compare results 
we have to renormalise the bulk cosmological constant, boundary cosmological constant 



and M to be consistent with [11]. In order to avoid doing this we note that the effect of 
such a renormalisation will to be to change the amplitude by a factor dependent only on p 
and q and so will not affect the ratio of correlation numbers. We therefore shall compare 



the ratio of correlation numbers. Computing the large w expansion of ( |A.1| ) we get, 

oo 

u> M (jt B ,lj) (1) =4^(-ir +1 (n-2)(n + 3)u;- 2 "- 1 (A.2) 

n=0 
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where the series converges if w > 1. We therefore find the torus 1-point correlation number 
for tachyon operators in the (2, 5) minimal string to be, 

, m An / im irn\ 2 „ / n\ „ / n\ , , n+i , , , , » . . , 
(^n) {1) = (sin — sin _) 2 r (l - -) T (--) (-1) — (n - 5)(n + 5)/i4 (A.3) 

when n is odd and zero if n is even. The 1-point correlation numbers on the torus for the 
(2, 2p + 1) models were computed in [11], giving, 



1 



(«„)"» = w> (1 > (-D^ ( ^] ' r ' 1 -^r> r 'i) (»+'-")(»+' 



sin 



r(i-sffiTO) 4p( P + i) 



2p+l / 1 V x 2p+l/ i Vn 



for n odd and zero when n is even. Using standard T-function identities it is then easy to 
show that the above equation is consistent with our computation (|A.3[ ). 

B. Liouville cylinder amplitude 

The full cylinder amplitudes were given in a usable form recently by [24|. We have obtained 
similar results independently and reproduce them here for the purpose of completeness. 
The cylinder amplitude between a (k, I) and (r, s) brane, Z{k, 1; &\\r, s; 02), is computed by 
the integral 

/•oo 

Z(k,l;a 1 \r,s;a 2 ) = drZ ghost (T)Z UouvilIe (r)Z matter (r), (B.l) 
J 

where r is the modular parameter of the cylinder in the closed string channel. The partition 
function for the ghosts is well known 

^ g ho S t(r) =r?(q) 2 , (B.2) 

where q = exp(— 2ttt) and r](q) denotes the Dedekind r/-function. The Liouville and matter 
contributions are given by 

7 _ S (k,l),(a,b)S(r,s),(a,b) ^ , . 

^matter - ff ~ ] X(o,6)W) K^- 6 ) 



(a,b)eE qtP 



(l,l),(a,6) 



roo 

liouville = / dP*t 1 (^)^ 2 (^)XP(q) J (B.4) 
J 

where E qtP is the Kac table for the minimal model, S is the relevant modular S'-matrix, 
and X( 0i 6)(q) and xp(q) are the characters for the minimal model primary field with Kac 
index (a, b) and the non-degenerate Liouville primary field respectively. The characters are 
given by, 

i 00 

X(a63 (q) = — — ^ (q( 2 P9 n +9°-P 6 ) 2 /4p9 _ ^pqn+qa+pbf/ipq^ ^ ^ ^ 



71= — OO 

,p 2 



Xp(q) =^r. (B.6) 
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We first concentrate on the ^matter- Recalling that the modular S-matrix for the (p, 
minimal model is, 



s, 



P 



(B.7) 



the explicit expression for Z ma _ ttei is 

■P-1 9-1 



-"matter 



V- E £ [(-D a(/+r+1)+6(fe+s+1) — ^ } sm( ^ } sin(7r2f) sin(7r2?) 



a=l 6=1 



sin(7rf)sin(7r^) 



»7(q) 



(2pqn+qa—pb)' 2 / ipq (2pqn+qa+pb) 2 /4pq 



(B. 



where we have used the symmetry of the Kac table to rewrite the limits on the a and b 
summation. Also, note that the quantity in the square brackets is symmetric in a and b 
and zero if a = or b = 0. This allows us to rewrite it as, 



j p— 1 9—1 
v ^ y o=l 6=-(g-l) 



(B.9) 



^^(cos-)^^^^)^!^-)^!^-)-^ £ q^ + *) 2 /^' 



where t = qa + pb. Substituting these expressions into ( B.2 ) gives, 
1 



Z 



°° dP 47f2 cos ( 27ro " ljP ) cos(2trj 2 P) 
sinh(2vr6P) sinh(^) 



/2pq 

['IX'fj 'TT'fj 'IX'fj IT'fj 'lY'fj ivfj 

sin( — ) sin( — )Uk-i(cos — )C/ r _i(cos — )[/;_i(cos — ){7 s _i(cos — ) 
v a v v a a 



o=l 6=-(g-l) 



dr ^ q (2p (i ,n+t) 2 /4p<J+ J P 2 



(B.10) 



Focusing our attention on the integral over the moduli, we compute, 

poo °° -i °° r -i 

dr £ q ( - 2pgn+t ^/ Apq+p2 = - ^ 

<n — — no in — — rv 



(2pgn + t) 2 /Apq + P 2 



1 



sinh( 



2-kP\ 



vi ■ 



2irP 



^\ZpqP cos — - cosh ■ , 



which gives, 



2tt 



2 /-oo dp cos(27rcT 1 P)cos(27rCT 2 P)sinh(^g) 2vriP 

Fk,l,r,s( — ;=) 

Jpq 



VlpqJo p sinh(27r6P)sinh( 2 ^; 



(B.ll) 



(B.12) 
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where F ki i^ s (z) is given by. 

F k,l,r,s( z ) = 



E E E E 

77=|jfc-r|+l,2 p=\l-s\+l,2 a=l 6=-(g-l) 



. 7Tt . Tit 

sm( — ) sm( — j x 



[/«-i(cos — ){J p _ 1 (cos — ) =7- 

P 9 cos ?l - cos z 



pq 



(B.13) 



where X p (k,r) = min(fc + r — 1, 2p — 1 — k — r). To compute ( |B.12|) we first note that 
the integrand is symmetric in P and so we can extend the integral to the whole real line. 
Without loss of generality we assume o\ > a 2 , and break up the factor of cos(27r(JiP) in 
( B,12| ) into its exponentials. For the integral containing, exp(27iwiP) we can close the 
integral along the real line by a semi-circle to produce the contour C+ in the upper-half 
plane, for the other term we use a contour, C_ in the lower-half plane. However, the 
integral using the contour in the lower-half plane, by a change of variables from P to —P 
is equal to the integral using the contour in the upper half-plane, we therefore get, 



Z(k,l;ai\r, s;a 2 ) = A c 



dP 



G(P,ai,cr 2 )Fk ) i,r > i 



,2vriF\ 



pq 



(B.14) 



The poles of F Khr>s {^) occur at P £ S F = { ^= : n £ Z, n ^ mod p, n / mod q}. 

hnP - 



The residue of F fc>1>P>s (2s£) at P G S F is, 



Res(F fc 



P 



in 



pq 



Jpq ,Kn ,irn. TT irn 
-2-± — r sm( — ) sm( — )Uf.-i{cos — )(/ r -i(cos — ) x 



2m 



P 



V 



P 



Trn 7rn 1 

Ui-i(cos — )C/ s _i(cos — ) . nn 
9 Q sm f? 



The residues of the G(P) are, 
27riRes \g(P):1 



27rfRes 
27uRes 
iRcs 



G(P);P 
G{P)-P 
G{P)-P 





< 




= ie 


— %<J\ 


inp ' 




n 




/pq. 




inq ' 




n 


2> 


fpq. 




inpq " 




-?-<- 

npq 


2> 


/pq. 





pq 



(B.15) 



(B.16) 



lfe T n (x 2 ) 



sin(7rn/g) 
sm(irpn/q) 
sm^nn/p) 
n s'm(irqn/p) 

) nip+q+1) e-^^T qn (x 2 ) 



Finally, a useful identity for computing the sums found in the computation of the 
cylinder amplitudes is, 



E? 7 



t=l 



cosh fit cos ( 



2a — ln(2(cosh 2 a + cosh 2 fi — 2j cosh a cosh fi cos <f> — sin 2 </>)) (B.17) 
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C. Loop equations 

To compute the resolvent Wx we make the following change of variable in (4.1), 

:i-a) 



S = Mi + M 2 + 1 + c, 



X 



, , (1 + a) „ , .a 



The partition function in these new variables takes the form, 



Z 



[dSdX] exp 



—TV (X 2 S + aXS 2 + V(S)) 
9 



where, 



V(S) = ^(1 + 3a 2 )S 3 - |s 2 + i(3c - 1)(1 + c)S 



(C.l) 



(C.2) 



(C.3) 



This explains why we chose the parameterisation in which Wm 1 +m 2 1S awkwardly at a = oo; 
we can compute it by computing Ws instead. As one would expect this matrix model is 
equivalent to an O(l) model, as demonstrated by making the change of variable, X — > 
(Z — aS)/2 which results in, 



J [dSdZ] 



exp 



■— Tr (Z 2 S/A + U(S)) 
9 



(C.4) 



where all dependence on a is absorbed into the new potential U(S). This shows that any 
correlation function that is odd in Z = 2X + aS is zero, which will be of use later. 

The resolvents can be obtained by generating an appropriate set of loop equations. We 
begin by obtaining the resolvents Wx and Ws- Consider the following change of variables 
in the matrix model, 



X ^X + e 



1 



z-S 



this gives the loop equation, 



W x s (z) 



a az 



W s (z) 



(C.5) 



(C.6) 



Now consider, 



Transformation 
1 



X -+ X + e 
S -> S + e 
X ^ X + e 

together with 



l l | l l 

x-X z-S ~f z-S x-X 

l l | l l 

x-X z-S ~r z-S x-X 
T^Sx^X -z-S + 

and ( p.6|) gives, 
1 



(C.7) 
(C.8) 
(C.9) 



9 



^x( z i x ) = ^s(z, x)VL S x(z, x ) + j-p ( w sx;s(z, x\ z) - aW S X;x(z, x\ x)) (CIO) 
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where we have introduced the following functions, 

n s (z, x) = W s (z) - aWx(x) - -(x 2 - a 2 z 2 + V'(z)) 

9 

n sx = w sx (z, x ) + -(x-—) (en) 

9 2 

n x = W?(z, x) + f + i(x - f ){x 2 - a 2 z 2 + V'(z)) 
2 g 2 

where P = Q - a 2 S + ax - ^ and Q = V'(z) - V'(S). 
Finally, the third transformation ( |C.8| ) gives, 

n S x(z,x)n S x(-z,x) = - (w s (z) + W s (-z) + aWx(x) + -(x 2 -^)) (C.12) 
- jpWsx;Sx{z,x;-z,x) 
Eliminating £lsx(z,x) between this and ( CIO ) gives, 



-n^(z,a;)n$(-z,x) = ns(«,x)ns(-« 1 x) ^) + if s H + *(x) + -(x 
9 V 54 

~ ]u2 (^5(2, 2;)^5X;5x(^, ar; ~z, x) + £l x (z, x)W S X;s{-z, x; - z) (C.13) 

-atl x (z, x)W S X;x(-z,x;x) + Q x (-z,x)(w S X;s(z, x; z) - aWsx;x{z,x;x)X\ 

By substituting in the large N expansion for all resolvents appearing in this expression 
we can obtain a recursive definition for the genus h resolvents. The genus zero resolvents 
satisfy: 

-n x (0) (z,x)n x {0) (-z,x) = (c.i4) 
9 

Qs(z,x)^Q s (-z,x)^ ( Ws i0 \z) + Ws^(-z) + aW x ^{x) + -{x 2 - 

v 54 

This is an important equation and all subsequent loop equations will be related to this 
through derivative-like operations such as the loop insertion operator. It is important to 
note that the LHS is polynomial in z and so if we expand the RHS in x about any point, 
the resulting Laurent expansion must have coefficients that are equal to a polynomial 
function of z. A convenient point to choose for this expansion is infinity as the definition 
of the resolvent then coincides with the Laurent expansion, whose coefficients are the yet 
to be determined quantities, (TrX n ). Thus we produce a number of equations containing 
Ws^(z), 

dx —n s ( z ,xf> n s (-z,x)^x (c.i5) 



C o 



2nix 



Ws {0) (z) + W s {0) (-z) + qW x (0) (x) + -(x 2 - ^) ) ,,J : ) 

9 4 
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where Coo is a contour around infinity, p n (z) is a polynomial in z with a finite number of 
unknown coefficients, n is an integer and the integral merely picks out the nth coefficient 
of the expansion. From the large x behaviour of Wx(x) ~ - it is easy to show that for 
n > 6 the LHS is zero and for n > 4 the LHS yields an expression with no dependence on 
W s {0) . The non-trivial cases occur for n = and n = 2 which generate enough equations 
to solve for Ws^°\z). The large z expansion of equations obtained for n < give relations 
among the quantities (Tr5 n ). The explicit equation obtained for Ws^°\z) is, 

Ws^(zf + ^(^ + V'(-z)-2V'(z))w s ^(z) 2 (C.16) 

+ ^(V(-z)(3z 2 a 2 - 4V'(z)) - 3z 2 a 2 V(z) + W (zf + P 2 {z)) W s {0) (z) 
1 



+ -^(Po(z) + (z 2 a 2 - V'(z))P 2 (z)) = 



4.9 



where, 



p t ( z ) = - Pt (z) + (Coefficient of (Ws (0) (z))° in pT5| for n = i) (C.17) 



The unknown constants in p n (z) are not all independent and they may be found in terms 



of the constants (TrS n ) by expanding the RHS of ( C.15 ) about z = oo and equating p n to 
the polynomial part of the Laurent expansion. The requirement that the singular part of 
the expansion vanishes again gives relations between the quantities {TrS n ). 



To find Wx^°\ we again consider ( |C.14 ), however we now know in terms of a 



finite number of unknown constants. On the LHS we have the function Cl x ^°\z, x) where 



P is a polynomial in z and S as defined in ( |C.ll ), whose highest power for both z and S 



is d — 1 and so it contains d unknown function of the form W^™ (x) for < n < d. 

We can generate a system of equations for these resolvents by again expanding both 
sides of the loop equation about infinity but this time in terms of z. This results in, 



dz 

z n+l 



-{l x i0 \z,x)n x i0 \-z,x) - (C.18) 
9 



n s (z,x)^n s (-z,x)^ (Ws^(z) + W s (0 \-z) + aW x ^{x) + -{x 2 - ^±.) 







Again we may vary n to generate different equations, however we now generate non- 
trivial equations for all < n < d. Which is enough to solve for all unknown resolvent 
functions appearing in 0, x ^°\z,x). However, unlike Ws^ for which we obtained an equa- 
tion for general V(x), we were unable to do this for the equations for W x "(x). We also 
note that for the equation giving Wx , the highest power of Wjf' ' appearing depends 
on the order of V{x). 

The approach thus outlined gives an algorithmic solution to this matrix model with 
an arbitrary potential and results in algebraic equations for all resolvents. 

We now apply the above procedure to the case when the potential is ( |4.12| ). There 
are two unknown constants in Pq(z) and P2(z) corresponding to (TrS) and (TrS* 2 ). Their 
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values are determined by the requirement that ( |C.16|) is a genus zero curve - which, in this 
case, is equivalent to the one-cut assumption. 

For Wx we can generate two equations from (C.18), for n = and n = 1, containing 
WOf' ' and W x ^ . The equation for Wj' ' resulting from the elimination of has a 

number of interesting properties. Generically it is quartic in Wx ■, however it reduces to 
a cubic for a = ±1 and it reduces to a quadratic of (Wx^ + polynomial in x) 2 for a = 0. 
Because of this the scaling limit for generic a does not hold for these particular values of 
a and must be taken separately. 

As was discussed before, the resolvents corresponding to the spin + and — Cardy states 
are Wm^ and Wm 2 respectively, which may be computed directly from (4J) and are 
clearly identical. However a faster method, since we already have Wx > is to set a = — 1 
implying X = Ah + |(1 + c). 

It is easy to generate loop equations for a general 2-loop amplitudes of the form Ws-H 
and Wx-h where H is a string of S and X matrices of length jVh- We can write H as, 
H = n^=i Xni where Xn is a matrix defined by, 



and Ix and 1$ are disjoint indexing sets. For such a product, the notation Huj\ is denned 



as , 

3 



H = Hxi (C20) 



The loop equations for Ws-h, can be generated by the following changes of variables. 

X -> X + e —^—TiH(S, X) (C.21) 

z — S 

X + 1 ' 1 +h.c)ltH(S,X) (C.24) 
\z\ — b x — X Z2 — & J 

The loop equations generated by X — > X + eF(X, S)TrH, may by obtained from the loop 
equations generated by X — > X + eF(X, S) by the following procedure, 

W Gi(X,S) -> W Gi(X,S);H 
W Gi(X,S) W G 2 (X,S) -> W Gi(X,S) W G 2 (X,S);H + W G2 (x,S) W Gi(X,S);H 

-> + E Wfyx.o^.S)* (*) (C.25) 



where Gi(X,S) is any function of X and S 1 . An equivalent set of rules applies for loops 
equations generated by S — > S + eF(X, S)TicH. The resulting loop equations can then be 
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solved by using the same method as used to solve the one- loop amplitude equations. For 
example, the two loop amplitude Ws-s( z ; u ) ma y be obtained by choosing H = (u — S)^ 1 
in (|CT2l| ), (|CT23D and (jOgg ). Applying the rules in (|CT^5| ) we get, 



and 



^5X(^, x)tt S ;S(z, X) U) + <9 U 



W sx (u)-Wsx(z) 



u — z 



(C.26) 



flsx(z,x)Q S X;s(-z,x;u) + r2 S x ; s(^,a;,'w)r2s'x(-^,a:) 

1 



.9 



(W S ;s(z; u) + Ws^-z; «) + aW X ; S {x; u)) (C.27) 



By eliminating x, u) we obtain an equation relating Ws ; s(z; it), Ws ; ,g(— z; u) and 

Wx-s{ x 'i u )- If we arrange the resulting equation so that the LHS is polynomial in z, we may 
expand the RHS in x about x = oo to generate a number of equations for Ws- t s(z', u) and 
Ws-s(—z; u). This is exactly the same procedure we used to compute the disc amplitudes 
in which the resulting expression contained a number of unknown constants, which for 
the equal potential case were (TrS) and (TVS' 2 ). However here we get an expression for 
Ws ; s(z', u) in terms of unknown functions of the form W s ; s{u) and W s -s( u )- These 
unknown functions may be found by computing the expansion of Ws- t s at z = oo and, 
by the symmetry of Ws ; 5, equating it to the expansion at u = oo. The expressions 
for W s -s{z) and W s2 ' -s( z ) then depend on a number of unknown constants of the form, 
(TrS' m TrS' n ) c , not all of which are independent, as can be seen by considering the large z 
expansion of the W s -s(z) and W s<2 -s{z). In the end the independent unknown quantities 
correspond to, (TrS'TrS') c and (TrS'TrS' 2 ) c which can be fixed by requiring that W s - S {z) 
has no singularities besides those at the branch points of Ws- The loop equations then 
give Wx-s i R terms of Ws-S- Given that the continuum limit for generic values of a does 
not give any more information we choose to set to a = —1 without loss of generality. 
Computing the scaling form of Ws-x gives ( |4. 18 ) and the other two loop amplitudes are 
found in a similar manner 

In order to compute the 1/N corrections it is clear from ( C.13 ) that Wsx-,S, Wsx-,x 



and Wsx-,sx need to be calculated. The method for calculating Wsx-,s was given in the 
preceding section. The calculation of Wsx-,x and Wsx-,sx, follows a similar procedure 
however the resulting loop equations contain new amplitudes of the form Wsx..SX- The 
calculation of these new quantities does not pose a great difficulty with the appropriate 
change of variables giving loop equations very similar to the ones presented previously. 

Once these quantities have been calculated and substituted in ( C.13| ), we may again 
expand both sides around x = oo and then z = oo to find an expression for Ws(z)^ 
and Wx(x)^\ Again this procedure introduces a new unknown constant corresponding 
to the -^r correction to (TrS) which may be determined by requiring that W s {z) {1) does 
not possess any poles besides those at the branch points of Ws(z)^°\ The result of the 
calculation is ( |4.22 ) 
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